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ABSTRACT 

We re-derive the equation for vertical hydrodynamical equilibrium of 
stationary thin and slim accretion disks in the Kerr spacetime. All the previous 
derivations have been unsatisfactory yielding unphysical singularities. Our 
equation is non-singular, more general, and simpler than those previously 
derived. 

Subject headings: accretion, accretion disks - - black hole physics - 
hydrodynamics — relativity 
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1. PHYSICAL MEANING OF THE VERTICAL EQUATION 



In this paper, we derive, in Kerr geometry, the vertical equation for stationary thin and slim 
accretion disks. The equation governs the balance between the vertical pressure gradient, 
vertical acceleration, and vertical components of all the inertial forces. Vertical direction 
means here the direction orthogonal to the plane of the disk, for example, the z direction 
in cylindrical coordinates, or the 6 direction in spherical coordinates. We assume that the 
disks are stationary and axially symmetric and that the vertical thickness of the disk H is 
much smaller than the corresponding cylindrical radius r. 

We derive the equation first in Newtonian theory in order to stress some points that are 
often overlooked, but the new result presented in this paper is the fully self-consistent 
derivation of the vertical equation in general relativity, in particular for the Kerr spacetime. 

In the reference frame comoving with the matter, the vertical pressure gradient force 
balances the combined vertical components of all the inertial forces: gravity, centrifugal 
force, the Euler force due to the vertical acceleration, and (in the case of a rotating black 
hole) the general relativistic Lense-Thirring force (dragging of inertial frames). The Coriolis 
force always vanishes in the comoving frame of reference and does not need to be considered. 

Previously, the vertical equation has been derived by Novikov & Thorne (1973; NT), Riffert 
& Herold (1995; RH) and Abramowicz, Chen, Granath k Lasota (1996; ACGL). None of 
these derivations was satisfactory. They all made incorrect assumptions that introduced 
unphysical singularity into the final form of the vertical equation. Formulae derived by 
NT and RH are singular at the location of the circular photon orbit, while the formula 
derived by ACGL is singular at the horizon. NT based their derivation on a discussion of 
the gravitational tidal forces in the vertical cylindrical direction. This set a standard for 
several further derivations, but (as pointed out by Lasota & Abramowicz, 1996) was not a 
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fortunate choice, because a mathematically correct discussion of the vertical tidal forces at 
the horizon (not yet given) would be very complicated and thus impractical. On the other 
hand, a direct approach based on the relativistic Euler equation is straightforward and 
simple. Here, we derive the equation directly from the relativistic Euler equation and make 
no additional simplifying assumptions. Our derivation is based on a general discussion 
given by Lasota & Abramowicz (1996). 



Physical quantities X(r, z) = X(r, —z) that are symmetric with respect to the plane of the 
disk, z — 0, are expanded around the plane according to 



Physical quantities antisymmetric with respect to the plane X(r,z) = —X(r, —z) have 
expansion 



Putting these expansions into the partial differential equations that describe stationary 
accretion disks, one converts the partial differential equations into set of N = oo equations 
with N = oo unknowns. Truncating the expansions at N — 0, thus keeping only zeroth 
order terms in (z/r), would not give a closed set: in this case there are more unknown 
functions than the equations. Truncating at N — 1 gives a closed set, with the number of 
unknown functions equal to the number of equations. The same is true for any N > 1. 
Only terms of the order (z/r)° and (z/r) 2 appear in the N — 1 equations. Therefore the 
approximation based on truncating all expansions at N — 1 offers two advanteges: (1) it is 
mathematically well defined and consistent, (2) it involves the smallest number of equations 



2. 



THE METHOD: SLIM ACCRETION DISKS 



X(r,z)= E X*(r)(-) 
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to be solved. The precise mathematical meaning of the "slim accretion disks" approach 
(Abramowicz et al. , 1988) is just that: one should keep all the terms of the order (z/r) 2 
and (z/r)° in all the equations. This should be contrasted with the "thin accretion disks" 
approach, where terms of the order of (z/r) 2 are kept in some of the equations, but in others 
they are rejected according to the tradition established in the influential paper by Shakura 
& Sunyaev (1973). It is now well understood that this traditional approach often causes 
serious physical inconsistencies which make the thin disk models inadequate for describing 
several astrophysically interesting types of accretion disks. 

Here, we derive one particular equation — the vertical equation — which belongs to the 
set of the slim accretion disk equations. All the other equations have been derived in the 
Kerr geometry by Lasota (1994) and Abramowicz et al. (1996). The vertical equation 
is itself of the order (z/r) 2 . We keep all the terms of zeroth and second orders in our 
derivation. Thus, our final "slim disk" equation contains all the (z/r) 2 terms that follow 
from a mathematically consistent procedure. However, we also discuss here the "thin disk" 
form of this equation in which some of the (z/r) 2 terms have been rejected in a way which 
is not mathematically consistent, but follows the traditional thin disk approach. 



3. NEWTONIAN DERIVATION 



It will be helpful to derive the vertical equation in Newtonian theory in both cylindrical and 
spherical coordinates. Of course, the final result does not depend on the coordinates used. 
However, there are some interesting differences between the equations written in cylindrical 
and spherical coordinates. When we use cylindrical coordinates, there is no centrifugal 
force in the "vertical" direction, while when we use the spherical coordinates, there is no 
gravitational force in the "vertical" direction. As we will see, it is exactly this property that 
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makes the spherical coordinates much better adapted for descring the flow near the black 
hole horizon. All previous derivations have been carried out in cylindrical coordinates, and 
we feel that it would be proper if we explain our derivation in these coordinates first. 

3.1. Cylindrical coordinates 

Let H = z(r) describe the location of the disk surface in cylindrical coordinates (r,z,<f>). 
According to our basic assumption that terms up to the second order should be kept, the 
pressure P(r,z), in terms of (z/r), has the expansion, 




where the small dimensionless parameter (z/r) gives the distance from the plane of the 
disk. The minus sign before (z/r) 2 P 2 (r) has been chosen to make P2(r) > 0. Po( r ) is the 
pressure at the plane of the disk. At the disk surface, P(r, z) = 0, and therefore 




Note that for z ~ H higher order terms in expansion (1) are not numerically small compared 
with the lower order ones. This is a weak point of the method but we have checked that 
the numerical inaccuracies introduced by this are of the same order as those of other 
approximate methods. 

From equations (1) and (2) it follows that, 




We shall denote by po( r ) the density at the plane of the disk. We do not assume anything 
about the vertical distribution of the density of matter p(r, z) as all derivations based on the 
"vertical integration" do. In particular, we do not assume that the density vanishes at the 
disk surface, or that it is in accordance with a polytropic equation of state, P = Kp 1+1 ^ n , 
as assumed by Hoshi (1977). Hoshi's formula reduces to our equation (3) for n — 0. 
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We shall write for the radial and cylindrical components of velocity, 



v r = v (r) + O 2 - 



(4) 
(5) 



For stationary flows, the vertical equation has the general form, 



1 (w) + (£' 

p \OZ I \oz 



dv z dv z 
or oz 



0. 



(6) 



Here $ = Q(R) is the gravitational potential that depends only on the spherical radius 
R = \fr 2 + z 2 , (dR/dz) = z/R. The Keplerian orbital velocity V K (R) in the potential 
is given by, 



Taking these relations into account, we write equation (7) in the form, 



(7) 



-2-+ - 

p \r J 



2 r 



T ri dux o 

V K + rv — h u ± - uiv Q 

dr 



+ o 6 [- 

,r 



0. 



(8) 



Thus, by demanding that the expression in the curly brackets vanishes, we find the most 
general version of the vertical equation, 



2* + (*V 

Po V r J 



V K + rv -— +u 1 - uiv 
dr 



0. 



(9) 



In this most general version, an additional unknown function Wi(r) appears, that is not 
present in the standard thin or slim set of equations. The additional equation which closes 
the set (not needed in the standard treatment where u± — is assumed ad hoc) is provided 
by the surface boundary condition which says that there should be no component of the 
total velocity orthogonal to the surface of the disk. In cylindrical coordinates this means 
that 

(10) 



Ui 

v 



dlnr 
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^From the above equation it follows that U\ and v must be of the same order with respect 
to (H/r) and that U\ = 0, is only a good approximation if H varies very slowly with r. 

Therefore, in the most general situation, the vertical equation consists of two equations: 
equation (9) which governs the vertical balance of forces and acceleration, and equation 
(10) which demands that in a stationary disk there is no component of velocity orthogonal 
to the surface. 

From equation (9) it follows that in general -Po/Po ~ Cf ~ (H/r) 2 with Cs being the sound 
speed. In addition, it is obviously V£ ~ (H/r)°. One can make simplifications in equation 
(9) by considering the (H/r) order of each term. In the standard thin disk model the radial 
flow is very sub-sonic, and therefore Vq ~ U\ <C Cg- Therefore, the terms rvo(dui/dr), u\, 
and U\Vq are of order higher than (H/r) 2 and may be dropped. The standard thin disk 
version of the vertical equation therefore has the form, 



In the transonic part of the slim disks the terms containing u\ are exactly of the (H/r) 2 
order and may also be neglected. From consideration of the vertical equilibrium alone it 
seems possible that the thin disk version of vertical equation could be also written in the 
form (11), without the terms containing u±. However, if these terms have been dropped, 
then in doing calculations with the full set of equations, one should perform a consistency 
check, and calculate the ratios (uiV )/Cg and (uiUi)/Cg using eqn. (10). If these ratios are 
~ 1, one must use the general form of the vertical equation, provided by equations (10) and 
(9). This is the case close to the horizon of the black hole where v ~ u\ ~ c 3> Cs- We 
shall return to this point later. 




(11) 



3.2. Spherical coordinates 



In spherical coordinates (R, 9, <fi) the small parameter that characterizes the distance from 
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the plane of the disk should now be cos 6, and the equation for the surface of the disk is 
given by 9 = Q(R). Obviously, in this case H/r xs cos©. According to this, a formula for 
the pressure, similar to equation (3) should be written as, 



P(R,9) = P (R) 



( cos 2 9 



IdP 2^cos^ 
and -— = -^-P (R). (12) 



v cos 2 

For the radial, azimuthal and polar components of the velocity we again write, 

v R (R,9) = v (R) + O 2 (cos 9), (13) 

V<P (R,9) = v^(R) + 2 (cos9), v e (R,9) = cos 9 Ul (R) + O s (cos 9). (14) 

Note, that does not have the physical dimension of velocity. This is because, in 
spherical coordinates, = (d<f)/dt) = Q, with Q being the angular velocity, and 

= g^v^ = R 2 sin 2 6VL = r 2 Q = £, with C being the angular momentum per unit mass. 
Similarly, v e has the physical dimensions of angular velocity, and both vq and U\ have 
the dimensions of angular momentum per unit mass. Of course, v R , Vr and v have the 
dimensions of velocity. One should bear this in mind when checking the physical dimensions 
of our formulae in spherical coordinates, or when comparing our formulae with those given 
e.g. by Tassoul (1978). Tassoul's azimuthal velocity "t^", which we denote here (to avoid 
confusion) by V^), is given by = i?sin 9v^ = v<p/Rsm 9, and has the physical dimensions 
of velocity. 

In spherical coordinates, the condition that there should be no component of velocity 
orthogonal to the surface of the disk, similar to that given by equation (10), takes the form, 

1 Ul _ din |f - 8] 



Rv d\nR 



(15) 



The general form of the vertical equilibrium is now given by a formula similar to (7), but 
written in spherical coordinates, 

1 dP „ divo/R) vgd(v e /R) v R v e vl . . 

- p m +Rv " A f 1 + i A W 1 + ^ + W cote = - (16) 
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From these formulae, using the same procedure as before, we derive the general form of 
vertical equation in spherical coordinates, similar to its cylindrical version (9) 



Note a significant difference between the two versions of the same vertical equation, written 
in cylindrical coordinates (9), and in spherical coordinates (17). While in (9) the Keplerian 
velocity V£ appears, in (17) we have the rotational velocity VZy The difference is due to 
the fact that in the z direction there is gravity force (~ V^/R) but no centrifugal force, 
while in the 9 direction there is centrifugal force (~ V^/R) but no gravity. Does this 
matter? One may argue that it does not, because 



and therefore, with accuracy to the higher order terms, these two versions of the vertical 
equation are identical (for the standard thin disk V£ = V^). 

But actually it does matter a lot. Firstly, in both Newtonian and general relativistic 
theories, the stationary accretion flows calculated in realistic 2-D and 3-D simulations 
are much more resembling quasi spherical flows (0 ~ constant) than quasi horizontal 
flows (H pa constant). This has been noticed by several authors, e.g. by Papaloizou & 
Szuszkiewicz (1994), and Narayan & Yi (1995). Thus, while u± — may be a reasonable 
approximation in spherical coordinates, it is not so in the cylindrical ones. Secondly, 
anticipating the relativistic discussion in the next section, the Keplerian velocity Vk(R) 
is singular at the location of the circular photon orbit, while the azimuthal velocity is 
non-singular everywhere, including the horizon. Thirdly, we know exactly what happens 
on the horizon in spherical coordinates. A general theorem (Carter, 1973) ensures that 
in spherical coordinates, U\(Rh) = 0. The use of cylindrical coordinates to describe the 
flow brings artificial, and technically complicated, difficulties near the horizon. These are 
generated only by a bad choice of coordinates (as, e.g. in Riffert & Herold, 1995). 




(17) 




~ — = cos B, 



(18) 



r 
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4. RELATIVISTIC DERIVATION 



We shall now derive the vertical equation in Kerr geometry, using the same method as in 
the previous Section. It is convenient to work at the beginning with a general metric, 

ds 2 = g tt dt 2 + 1g H dtd$ + g w d4> 2 + g RR dR 2 + g ee d6 2 , (19) 

and specify the Kerr metric functions g ik = gik(R, 0) at the end of the calculations. 
The stress-energy tensor for a perfect fluid has the form, 

T\ = pu*u k - 5\P, (20) 

where u % = dx % /ds is the four- velocity of the fluid, and p = e + P, where e is the energy 
density (we use units in which c—l — G). The four- velocity is a unit vector, which means 
that 

1 = g ik UiU k = g a u t u t + 2g t<t> u t u 4> + g^u^ + g RR u R u R + g 09 u e u e . (21) 
The inverse metric g tk obeys g tk gjk = 

We derive the vertical equation taking the j = 9 component of the general equation 
h % jV k T\ = 0. Here is the covariant derivative operator, and h l ■ = S l j — u' l Uj is the 
projection tensor. This gives, at the end of simple standard calculations, 

ldP_ k du^_ ri k (22) 
P d9~ U dx* 1 6kU%U ' {22) 

where T l ek is the Christoffell symbol that should be computed from the metric components 
and their first derivatives. This relativistic equation provides the replacement for the 
Newtonian equation (16). As in the Newtonian case, the leading order of this equation is 
(cos^) 1 . Thus, in calculations that follow, we shall keep only these terms, and neglect the 
O s (cos8) and higher order terms. 
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The —Y l gk UiU k term should be computed with the help of equation (21). The result in Kerr 
geometry is 



4- n^V* 4- rj*^r* 4- rj^T* — ?n t(/> r R 



-u t u t 



g tt Y\ e + g*l* - g tt T 



^R 

Re- 



(23) 



In deriving this equation we took into account that uqUq ~ cos 2 6* should be dropped, and 
that in Kerr geometry g RR T e eR + g e8 V R e = 0. The term u k due/dx k equals 



k du e 
u ——r = cos V 



dx k 



g RR V ^- - g ee Ul 
y u dR y 1 



(24) 



where U\ and Vq are defined in terms of the following expansions [cf. equations (13) and 
(14)], 

u R = V (R) + 2 (cos6), u e = cos eU^R) + £> 3 (cos#). (25) 
One also needs to use the relations u 9 = g m ua and u R = g RR UR. 

The terms containing pressure gradient should be treated the same way as in the Newtonian 
derivation. 



To the desired orders, the components of the inverse metric, and the Christoffel symbols 
which appear in (23) and (24) are in Kerr geometry, 



9 tt = 

g<t>4> — 

9" = 

-p4> _ 

L te — 

r 4> _ 

J- MO 



2a 2 (R+2M)+R 3 
RA ' 

R-2M 



t(j> _ 2aM 
9 RA ' 



RA 

1 

' R 2 ' 
2aM 



cos 9, 



r* - 
1 te — 

J- AO — 



_A_ 

R 2 ' 

•^costf, 



R 3 

1 + ^W T R e 



2a? M 
R 3 



COS 9, 



■jp cos6> 



By inserting these into equations (22 - 24), we obtain the final result, 



Po 



R) \R 2 



~dR 



0. 



(26) 



(27) 
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Here C 2 = C 2 — a 2 {£ — 1), where C = —u^ is the conserved angular momentum for a 
geodesic motion and that £ = u t is the conserved energy for such a motion. Realistic 
matter almost free-falls (moves along geodesies lines) when it crosses the horizon. Thus, 
~ const near the horizon. 

Formula (27) above gives the most general version of the vertical equation for stationary 
flow. On the horizon, U\ = 0, and £* is finite, Vb — > oo, but A — > 0, in such a way that VoA 
is a finite quantity. Thus our equation is non-singular on the horizon. 

As in Newtonian theory, equation (27) should be considered together with equation (15) 
which, obviously, has the same form in relativity. However, in the spirit of the "thin disk 
approach" , we suggest that in practical applications, involving stationary flow, it is always 
safe to drop the terms containing U\. As we have already explained, they cannot be 
significant anywhere in the flow: far from the black hole because the flow there is sub-sonic 
or transonic, and close to the hole because flow there is quasi spherical. Thus, we conclude, 
that the vertical equation could be taken in practical applications in the form, 



The Schwarzschild version of this equation was first derived by Lasota and Abramowicz 
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